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ABSTRACT 

To test the theoretical understanding that finding bright CO emission depends primarily on 
dust shielding, we investigate the relationship between CO emission (Ico) and the amount of 
dust (estimated from IR emission and expressed as “Ay”) across the Large Magellanic Cloud, 
the Small Magellanic Cloud, and the Milky Way. We show that at our common resolution of 
10 pc scales, Ico given a fixed line-of-sight Ay is similar across all three systems despite 
the difference in metallicity. We find some evidence for a secondary dependence of Ico on 
radiation field; in the LMC, Ico at a given Ay is smaller in regions of high Tdust, perhaps 
because of an increased photodissociating radiation field. We suggest a simple but useful pic¬ 
ture in which the CO-to-H 2 conversion factor (Xco) depends on two separable factors: (1) 
the distribution of gas column densities, which maps to an extinction distribution via a dust- 
to-gas ratio; and (2) the dependence of Ico on Ay. Assuming that the probability distribution 
function (PDF) of local Milky Way clouds is universal, this approach predicts a dependence 
of Xco on Z between Z~^ and Z~^ above about a third solar metallicity. Below this metal¬ 
licity, CO emerges from only the high column density parts of the cloud and so depends very 
sensitively on the adopted PDF and the H 2 /H 1 prescription. The PDF of low metallicity clouds 
is thus of considerable interest and the uncertainty associated with even an ideal prescription 
for Xco at very low metallicity will be large. 
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1 INTRODUCTION 

As the immediate reservoir for star formation, the molecular inter¬ 
stellar medium (ISM) plays a key role in the evolution of galax¬ 
ies. Unfortunately, the majority of cold molecular hydrogen (H 2 ) 
in typical clouds is invisible in emission due to the fact that the 
H 2 molecule has low mass and therefore requires high tempera¬ 
tures (the lowest level corresponds to E/k « 510 K) to excite 
its rotational transitions (Kennicutt & Evans 2012). As a result, as¬ 
tronomers employ a suite of more observationally accessible trac¬ 
ers of H 2 to study molecular clouds. Low J rotational transitions of 
CO represent the most accessible and commonly used such tracers, 
especially in external galaxies. The abundance, pervasiveness, and 
brightness of CO make it a useful tracer, but it is not perfectly co- 
extant with H 2 and the relationship between CO emission and H 2 
column density (so called “CO-to-H 2 conversion factor”; Xco) is 
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both expected and observed to vary systematically as local condi¬ 
tions change (see the review by Bolatto et al. 2013). Consequently, 
it is important to understand the physical origins of CO emission 
and their implications for the use of CO emission to trace H 2 . 

Coarsely, variations in Xco are twofold. The regions in a 
molecular cloud where CO and H 2 exist are not perfectly matched, 
with a layer of CO-poor H 2 extending beyond the region at which 
the dominant form of carbon changes from CO to ClI (e.g., 
Maloney & Black 1988; van Dishoeck & Black 1988). Theoreti¬ 
cally, the amount of dust shielding (Ay) between the CO-ClI tran¬ 
sition layer and the H 2 -H 1 transition layer is estimated to be almost 
constant (Wolfire et al. 2010). This is also seen in PDR calculations 
(Bell et al. 2006) and numerical simulations (Glover & Mac Low 
2011; Shetty et al. 2011), where Ico exhibits a clear dependence 
on Ay. In a low metallicity cloud the dust-to-gas ratio is also low, 
so that achieving some fixed Ay requires a much larger column of 
gas than at high metallicity. Therefore this intermediate region of 
H 2 without much associated CO becomes very large in terms of to- 
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tal gas content. As a result, one expects to find less CO emission per 
unit H 2 in clouds with a low metallicity. This gives rise to a dust-to¬ 
gas ratio, and thus metallicity, dependent term in Xco- Meanwhile, 
within region where CO is abundant, the line is usually optically 
thick. This leads the ratio of CO emission to gas mass in this region 
to depend on density, temperature, and potentially other dynami¬ 
cal factors (e.g., see Maloney & Black 1988; Downes & Solomon 
1998; Narayanan et al. 2012). 

In this paper, we focus on the first part of the problem, the 
metallicity-dependent term in Xco- We focus on the relationship 
between dust abundance along the line of sight and the brightness 
of CO emission (the “Ico-Av” relation) as a way to explore the 
physics of CO emission and its utility as a tracer of H 2 . Because 
theoretical models highlight the key role of dust shielding in set¬ 
ting the extent of bright CO emission, our hypothesis is that across 
diverse environments we will often find about the same amount of 
CO emission per unit dust shielding. 

Highly resolved (sub-pc) observations of individual nearby 
clouds have explored the Ico-Ay relationship in detail 
(Lombardi et al. 2006; Pineda et al. 2008, 2010). These provide 
strong observational support for the crucial role of Ay in deter¬ 
mining the amount CO emission. These high resolution studies 
also reveal distinct regimes in the relationship, such as an extinc¬ 
tion threshold below which CO emission is faint or absent, a linear 
rise of Ico at intermediate extinctions, and evidence for satura¬ 
tion at high extinctions. The relationship has not been explored as 
much outside the Milky Way because of the coarse physical resolu¬ 
tion in most CO and dust maps of other galaxies. However, study¬ 
ing a molecular complex in the Small Magellanic Cloud (SMC) at 
« 10 pc resolution, Leroy et al. (2009) did observe an Ico- Ay 
relationship that resembled that for a Milky Way cloud (their Fig¬ 
ure 7). 

Several new data sets make it possible to revisit the relation¬ 
ship between dust and CO emission in the Magellanic clouds over 
a much wider area. The Magellanic Mopra Assessment (MAGMA) 
(Wong et al. 2011) obtained high spatial resolution (~ 10 pc) CO 
J = 1 —0 data across most areas of bright CO emission in the 
LMC. A new APEX survey of the southwest part of the SMC (Ru¬ 
bio et al., in prep.) provides similar coverage in that galaxy; we also 
use previously published CO J = 2 —>■ 1 and J = 1 ^ 0 maps 
of the N83 complex in the SMC Wing (Bolatto et al. 2003). Key 
projects by Spitzer (Surveying the Agents of a Galaxy’s Evolution, 
SAGE Meixner et al. 2006) and Herschel (HERschel Inventory of 
The Agents of Galaxy Evolution, HERITAGE Meixner et al. 2010, 
2013) allow us to model infrared (IR) emission to estimate the line 
of sight extinction. This provides a handle on the total extinction (or 
dust column) through a part of the galaxy, which offers an imper¬ 
fect but observationally accessible analog to the physically crucial 
shielding of material from the interstellar radiation field (ISRF). 
Clearly, the relation between the total dust column (expressed by 
“Av”) and the degree of shielding towards an average CO molecule 
depends on geometry, but even in simulations Ay and real shield¬ 
ing appear closely related (e.g., see Glover & Mac Low 2011). 
Thus, we now have a handle on CO emission and dust column at 
« 10 pc resolution over a matched area for the two nearest (~ 
50 kpc for the LMC and 60 kpc for the SMC Keller & Wood 2006) 
star-forming low-metallicity galaxies (~ 1/2 Zq for the LMC and 
1/5 Zq for the SMC Westerlund 1997). 

Combining these data on the Magellanic Clouds with Milky 
Way data from Planck and the Dame et al. (2001) CO survey, we 
are able to ask how the CO intensity, Ico, at a given line of sight 
extinction. Ay, compares between the Milky Way, the LMC, and 


SMC. Do successive steps of a factor of « 2 in metallicity have 
a visible impact on Ico at a given Ay or is the amount of dust 
shielding alone the key parameter? We also search for secondary 
factors affecting Ico at fixed Ay, with the most obvious candidate 
being the interstellar radiation field, which is directly traced by the 
dust temperature. This might be expected to influence the amount 
of CO emission at a given line of sight extinction in two ways: first 
lowering the amount of CO emission by raising the number of dis¬ 
sociating photons and so requiring more dust shielding for the tran¬ 
sition from ClI to CO. Second, perhaps increasing the temperature 
of the CO and so increasing /co¬ 
in the second part of the paper, we explore the implications of 
a universal Ico-Ay relation for the CO-to-H 2 conversion factor, 
Aco- If Ico is largely set once Ay is known, then the distribu¬ 
tion of Ay becomes the key factor to predict CO emission. This, in 
turn, depends on the probability distribution function (PDF) of gas 
column densities and the gas-to-dust ratio. The PDF of individual 
clouds in the Milky Way has been the study of significant quantita¬ 
tive study in recent years (see Kainulainen et al. 2009, and follow¬ 
ing). We combine these results with our estimates of the Ico-Ay 
relation to make an empirically-driven estimate of how Aco de¬ 
pends on metallicity. 


2 DATA AND MODELLING 

We aim to compare CO emission to the line-of-sight extinction, 
estimated from IR emission, on the scale of individual clouds (~ 
10 pc). To do so, we assemble matched-resolution CO (Section 2.1) 
and infrared (IR) emission maps (Section 2.2) for the LMC, SMC, 
and Milky Way. We use the IR data to estimate the line of sight 
extinction (Section 2.3) and so estimate Ico as a function of Ay. 

2.1 CO Data 

2.1.1 LMC 

We use the second release of the Magellanic MOPRA Assess¬ 
ment (MAGMA Wong et al. 2011) survey^ to generate an inte¬ 
grated intensity (“moment 0”) map of CO emission from the LMC. 
MAGMA used the 22-m MOPRA telescope to observe the CO 
J = 1 —>■ 0 transition toward molecular clouds identified from 
NANTEN surveys (Fukui et al. 1999, 2008). MAGMA has an an¬ 
gular resolution of 45”, « 10 pc at the distance of the LMC. The 
mean rms brightness temperature of the MAGMA cube in a single 
channel (0.5 km s“^) is 0.3 K (Wong et al. 2011). We generate the 
integrated intensity map by directly integrating the cube along the 
whole velocity axis (180 km s“^ < Vlsr ^ 320 km s“^). 

Eor both theoretical and practical reasons, our analysis will 
treat Ay as the independent variable. Ay is expected to set the 
amount of CO emission and IR emission is detected at higher sig¬ 
nal to noise than CO throughout the Magellanic Clouds. Reflecting 
this, we work with integrated intensity derived from a broad veloc¬ 
ity window that will include any CO emission from the LMC. This 
directly integrated intensity map is not clipped or masked and so 
includes both positive and negative values. The advantage of this 
approach is that by averaging together many spectra at the same 
Ay we can recover a mean Ico that is too faint to be detected in 
an individual line of sight. A corollary of this approach is that we 
are sensitive to small zero point offsets in the data. Therefore, we 
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subtract a constant baseline from the MAGMA data cube pixel-by- 
pixel, with the value determined from the median intensity of each 
line of sight from the signal-free edge channels at the edge of the 
data cube. The median offset in the zero level is very small, corre¬ 
sponding to only 1.2 mK and only important because some of 
our analysis focuses on faint regions. In Appendix A1 we show the 
quantitative effects of varying our baseline treatment, which is very 
minor compared to other uncertainties in the analysis. 

2.i.2 SMC 

For the SMC, we use a new APEX survey of CO J = 2 —>• 1 
emission from the southwest region of the SMC data (P.I.: Ru¬ 
bio, Rubio et ah, in prep.). The angular resolution of the data is 
~ 28”, corresponding to « 8 pc at the distance of the SMC, and 
rms noise is 0.3 K in each 0.1 km s“^channel. These data target 
the southwest part of the SMC bar, which contains most of the on¬ 
going star formation and molecular gas in the SMC. We also use 
SEST observation of CO J = 2 —>■ 1 and J = 1 —>■ 0 emis¬ 
sion from the star forming complex N83 located in the wing of the 
SMC (Bolatto et al. 2003). These data have angular resolution of 
38"(J = 2^1) and 55"(J =1^0) and rms noise 0.1 K in a 
0.25 km s“^channel for both lines. 

As in the EMC, we directly integrate these data along the ve¬ 
locity axis to generate integrated CO intensity maps of the south¬ 
west SMC, picking a velocity range that covers the whole region 
and carrying out no other “masking.” The velocity ranges of inte¬ 
gration are 80 km s“^ ^ Vlsr ^ 160 km s“^ for the southwest 
SMC and 140 km s'^ ^ Vlsr ^ 187 km s"^ for N83. 


2.2 Infrared Maps 

We use observations of IR dust emission in the EMC and the SMC 
at four different wavelengths from the HERITAGE survey: 100 and 
160 nm images from the PACS instrument, and 250 and 350 /rm 
images from the SPIRE instrument. The angular resolutions of 
these maps are 7.7”(100 fim), 12”(160 /rm), 18”(250 /rm), and 
25”(350 /rm), so that the CO data set the limiting resolution for our 
analysis. 

Before proceeding to our analysis, we convolved the IR maps 
in the EMC to the 45"resolution of the MAGMA CO map. This is 
done by first using the convolution kernels of Aniano et al. (2011) 
to convolve HERITAGE maps to a common PSF (we used the 
Spitzer 160^m for comparison with Spitzer work) and then degrade 
them together to the 45" resolution of the MAGMA CO data. We 
place all EMC data on the same astrometric grid, which has pixel 
spacing of 15.6". 

We take a similar approach to match the resolutions and grids 
of IR maps in the SMC. Here we match the Herschel data to the 
28" resolution of the CO J = 2 —1 transition in the southwest 
region of the SMC, to the 38" resolution for the CO J = 2 —1 
transition in N83, and to the 55" for the J = 1 —>■ 0 transition in 
N83. Again all data are placed on a shared astrometric grid. 

The 1-(T noises of the EMC IR maps at our working resolution 
(45") are crioo ~ 2.3 MIy sr“^, (Tieo ~ 2.4 MIy sr“^, (T 250 ~ 
0.88 MJy sr~^, and 0-350 0.48 MJy sr“^. In our analysis we 

consider only regions with intensity at least three times these values 
in each band. We refer to the “EMC field” as the region that satisfies 
this S/N ratio cut in the IR maps, while the “MAGMA field” is 
defined as the region in the EMC field where the MAGMA survey 
mapped in CO. The MAGMA field is a subset of the EMC field. 


The noises in the SMC IR maps are similar but vary with reso¬ 
lution. We estimate the noise at each working resolution and again 
mask regions below 3-o- in the IR maps. As above, the “SMC field” 
is the region where IR maps have values greater than 3-o- at high 
resolution, while the “APEX field” and “SEST field” are the regions 
within the SMC field where APEX telescope and SEST telescope 
mapped CO emission. 

Note that a first order (linear) baseline has been already sub¬ 
tracted from the IR maps in the HERITAGE release to remove 
emission not associated with the Magellanic Clouds themselves 
from the maps (Meixner et al. 2013). We do not apply any fur¬ 
ther correction to account for Milky Way foreground dust emis¬ 
sion, though we do model a level of uncertainty in this subtraction 
by including a zero point uncertainty in our Monte-Carlo analysis 
(Section 2.3.3). 

2.3 Estimation of the “Av” Map 

We use the Herschel IR emission maps to estimate “Ay” along 
each line of sight through the EMC and the SMC. Here Ay refers 
to visual extinction, measured in magnitudes estimated from the 
optical depth at 160 pm, rieo- We calculate this by fitting a mod¬ 
ified blackbody to the measured IR intensities (Section 2.3.1) and 
then converting the dust optical depth to visual extinction following 
an empirical scaling derived from the Milky Way (Section 2.3.2). 
Because the dust emission is optically thin, “Ay" measured in this 
way will probe material along the whole line of sight and averaged 
over the substantial beam of the Herschel data. This is similar to ex¬ 
tinction mapping for nearby molecular clouds, which uses sources 
behind the cloud to create a large scale map. It differs from true ex¬ 
tinction mapping using internal sources in the Magellanic Clouds 
(e.g., stars), which will measure the extinction only part of the way 
through the galaxy. 

2.3.1 Modified Blackbody Fit 

We assume that the dust along a line of sight can be described as 
an optically thin (r <C 1) greybody at an equilibrium temperature 
Tdust, and wavelength dependence of dust optical depth is a power 
law with spectral index /3, i.e., tx oc (e.g., Draine & Eee 

1984). In this case the optical depth at 160 pm (neo) is given by 

^_ heo _ 

160 pm)’ 

where /leo is the observed 160 pm intensity and Bu(Tdust, A) is 
the intensity of a blackbody of temperature Tdust at wavelength A. 

Because we have dust emission intensities measured at 
four different wavelengths, in principle we can fit for the 
three unknowns /3, rieo and Tdust- Instead we adopt a fixed 
P = 1.5 as our fiducial value and fit for two unknowns 

T 160 and Tdust by minimizing the from observed IR inten¬ 
sities and the model IR intensities, taking the colour correc¬ 
tion^ for each filter into account. We fix /3 in order to mini¬ 
mize the uncertainties on the Tdust (and thus rieo) arising from 

^ See http://herschel.esac.esa.int/twiki/pub/ 

Public/PacsCalibrationWeb/cc_report_vl .pdf for 

PACS colour coiTection and http://herschel.esac.esa.int/ 
hcss-doc -11 . 0/load/spire_drg/html/ch05s07 . html for 

SPIRE colour correction. 
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the fact that P and Tdust are somewhat degenerate in space 
(Dupac et al. 2003). The adopted /3 above is a reasonable de¬ 
scription of the integrated SED of the LMC (Bernard et al. 2008; 
Gordon et al. 2010; Planck Collaboration et al. 2011a) and the 
SMC (Stanimirovic et al. 2000; Aguirre et al. 2003; Leroy et al. 
2007), and is intermediate in the range of plausible astrophysical 
values, 1.0 < /? < 2.0 (Draine & Lee 1984, Schlegel et al. 1998). 


2.3.2 Conversion to Av 

Our modified blackbody ht yields the optical depth at 160 ^m, rieo- 
We wish to phrase our analysis in terms of the line of sight V -band 
extinction, Av, which is the conventional unit expressing shielding 
in discussion of PDR regions. We translate rieo to Av via 


Av ~ 2200 r 160 • (2) 

We arrive at this conversion in several ways and take the dif¬ 
ferences among the estimates to indicate the uncertainty in the 
conversion. First, following Leroy et al. (2009), we infer rieo = 
2.44 X 10"^® (V(Hl) from DIRBE/FIRAS observations of the 
Galactic diffuse Hi (Boulanger et al. 1996). Then, we adopt the re¬ 
lation between colour excess of Solar Neighborhood stars and Hi 
column density from Ly a absorption, E{B — V)= W(H)/5.8x 
10^^ cm“^ (Bohlin et al. 1978). Taking a Galactic Rv = 3.1, we 
estimate Av ~ 2190 rieo- 

We compare this number to the Schlegel et al. (1998) conver¬ 
sion used to estimate Galactic extinction maps from IRAS data. Af¬ 
ter correction to a fiducial dust temperature, they find E{B — V) — 
(0.016 ± 0.004) Iu)o, where Iiqq is the 100 ^tm intensity after 
correction to a fixed 18.2 K temperature, and E{B — V) is red¬ 
dening measurement of background elliptical galaxies. Again tak¬ 
ing Rv = 3.1 and assuming /3 = 2 for Galactic dust®, the 
Schlegel et al. (1998) conversion corresponds to Av ~ 1939 rieo- 
More recently, Planck Collaboration et al. (2013a) find systemati¬ 
cally higher dust temperature and thus lower dust optical depth at 
high Galactic latitude using IRAS and Planck data. This is mainly 
due to the observed dust emission SED being flatter (/) « 1.59) 
than previous studies. This leads to a higher E{B — V) per dust 
optical depth (r), where the former is from SDSS quasars. Adopt¬ 
ing Rv = 3.1 and /3 from their study, the conversion is Av ~ 
3246 r160. 

This conversion can also be related directly to the mass ab¬ 
sorption coefficient of dust, k, following Hildebrand (1983). We 
consider neo = EduatKieo. In the Milky Way, Hildebrand (1983) 
suggests K 250 ~ 10 cm^ g“^, which we convert to 160 fim as¬ 
suming P = 2. Then taking the relation between E{B — V) and 
W(H) above and assuming a dust-to-gas mass ratio of ~ l-to-150 
(Draine et al. 2007) and Rv = 3.1, we arrive at Av ^ 1443 ri6o. 

To synthesize, calculations based on Milky Way studies sug¬ 
gest Av = (1400 — 3200) ri6o, with the mean conversion corre¬ 
sponding to Av « 2200 ri6o, which is the conversion we use in 
this paper. The fractional uncertainty above approximately matches 
that of Schlegel et al. (1998), and we fold it in to our Monte-Carlo 
analysis to account for the uncertainties in the conversion from neo 
to Av (See Section 2.3.3). 

There are several direct measurements of E{B—V) and Av in 


® In comparing our tiqq-Av conversion to other works, we use p assumed 
in each study rather than our fiducial /? = 1.5. 
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Figure 1. Comparison of our FIR-based Ay estimate to Ay inferred from 
stellar photometry (Zaritsky et al. 2004), UV spectroscopy (Welty et al. 
2012), and NIR colour excess (Dobashi et al. 2008) in the LMC. The log 
ratio of the former to the latter (logj^Q ylFlR/ylVV,Opt,NlR^ shown on 
the a:-axis, and the fraction of lines of sight is shown on the j/-axis. The his¬ 
togram for Zaritsky et al. (2004) is shown on top of the Welty et al. (2012) 
histogram, which is on top of the Dobashi et al. (2008) histogram. We ex¬ 
pect the Zaritsky et al. (2004) and Welty et al. (2012) Ay estimates to sam¬ 
ple only about half of the LMC, on average, and so scale these estimates by 
a factor of two in the plot. The dotted line indicates a ratio of unity. This 
lies close to the center of the sum of three histograms, indicating an overall 
agreement between our approach to estimate Ay from the FIR and direct 
measurements, though there is substantial scatter and systematic effects re¬ 
main. (A colour version of this figure is available in the online journal.) 


the LMC. Our neo map has complete coverage and high signal-to- 
noise compared to these maps, so they do not offer a viable replace¬ 
ment, but we use them to check our adopted conversion from rieo 
to Av (Equation 2). We compare to three data sets: the Av map es¬ 
timated by Zaritsky et al. (2004), which is based on photometry of 
individual “hot” stars; the compilation of spectroscopic measure¬ 
ments from Welty et al. (2012); and the Ay map inferred from the 
near-IR (NIR) colour excess method by Dobashi et al. (2008) us¬ 
ing 2MASS data. In all cases, we restrict the comparison to regions 
that have Av > 0.5 mag in our FIR-based Av map. Zaritsky et al. 
(2004) and Welty etal. (2012) report extinction to LMC sources 
that we expect, on average, to lie halfway through the galaxy. To 
compare to our map, which samples the whole column, we multi¬ 
ply those data by a factor of two to account for the difference in 
geometry. The Dobashi et al. (2008) map already accounts for the 
distribution of stars along the line of sight and attempts to report an 
integrated extinction. Figure 1 shows the resulting comparison as a 
histogram of the ratio between our map and these other estimates. 

Our FIR-based map yields lower Av, on average than stel¬ 
lar photometry or UV spectroscopy, after the factor of 2 scaling of 
the Optical and UV estimates. We find median « 

0.75 for the Zaritsky et al. (2004) data and median A^F^'/A^ « 
0.95 comparing to the Welty etal. (2012) data. In both cases we 
observe significant scatter, « 0.24 dex in the Zaritsky et al. (2004) 
case and « 0.33 dex comparing to Welty et al. (2012). Much of 
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the large scatter likely reflects the internal geometry of the LMC 
or, for Welty et al. (2012), the difference between our large beam 
and the single pencil beam sampled by their spectroscopy. On the 
other hand, the FIR-based Av yields somewhat higher Ay than 
the NIR color excess method. We find median « 1.6, 

again with large scatter (0.33 dex). Altogether, the sum of three his¬ 
tograms is approximately centered at a ratio of unity (dotted line). 
Given the stark difference in approaches to estimate the extinc¬ 
tion we view these comparisons as reasonable confirmation of our 
adopted Aw/neo and overall approach. We take them to confirm 
the « 50% systematic uncertainty in Av/neo discussed above. 

2.3.3 Estimating the Uncertainties in rieo and Av 

We are primarily interested in the average relation between Ay 
and Ico in the Magellanic Clouds. To estimate the uncertainties 
involved in the Ay portion of this relation, we adopt a Monte-Carlo 
approach. We use the LMC data as a basis to repeatedly re-calculate 
our Ay map while varying our assumptions across their plausible 
range. Each time that we do so, we add a new realization of the 
statistical noise to the data. We begin with the true LMC maps. We 
then add Gaussian noise with magnitude matched to the measured 
statistical noise to each map. For each new map, we also vary the 
zero point of the maps within the estimated uncertainty. Doing so, 
we generate 100 new maps that could be realistic observations if 
the observed LMC maps were indeed the true intensity on the sky. 
For each set of these noise-added maps, we generated 10 sets of 
Tduat and rieo maps following the minimization described in 
Section 2.3.1. For these maps instead of fixing /3 at 1.5, we fixed 
it randomly at a value within its plausible ranges, 1 < /3 < 2. In 
the end, we have 1000 maps with a spread that captures our true 
uncertainty regarding the derivation of Ay, except the conversion 
from ri6o to Ay. We take this into account by randomly taking 
a plausible conversion from rieo to Ay as determined in Section 
2.3.2, Av/ti 60 = (1400 — 32 00), resulting in 1000 Ay maps. 

Based on these calculations, we calculate the a 1-a fractional 
scaling uncertainty associated with Tduat, 10%, and with Ay, 45%. 
That is, the whole map is uncertain by this amount due to zero-point 
uncertainties and methodological decisions. This is a correlated er¬ 
ror that will adjust the entire map. The exercise and the scatter 
would be almost the same for the SMC, so we take these errors as 
representative of both galaxies. In the following figures discussing 
Ico-Ay relationships, we show this error estimate for Ay as the 
horizontal error bars in the bottom-right comers to represent the 
typical uncertainty in Ay. 


2.3.4 Limitations of Our Approach 

Our approach to estimate “Ay” has limitations, both due to our 
adopted approach and the use of IR emission to trace dust. We 
model a single population of isothermal dust along each line of 
sight. In reality, a mixture of temperatures and grain properties 
are present along each line of sight. This leads to biases in the to¬ 
tal dust column determination (e.g., see Schnee et al. 2005, 2006, 
2007, 2008 for detailed discussion) and could potentially affecting 
Ay /tiqq. These biases are somewhat alleviated by the inclusion 
of the long wavelength Herschel data and the very high (by extra- 
galactic standards) spatial resolution of Herschel at the LMC. Ulti¬ 
mately, they correspond to fundamental degeneracies in modelling 
the IR SED to derive a dust column. Resolution clearly represents 
another limitation; while 10 pc resolution is the best achievable 


Table 1. Comparison of Tduat and Ay in MAGMA field 



< Ay > 

(t(Av) 

'^dust ^ 



(mag) 

(mag) 

(K) 

(K) 

Spitzer 

0.911 

0.940 

21.8 

1.84 

Herschel 

0.637 

0.627 

24.1 

2.82 


outside the Milky Way, this is still very coarse compared to sub¬ 
structure observed in Milky Way clouds, so that measured “Ay” 
corresponds to something more like an average extinction across a 
Milky Way cloud than a value within a cloud. Finally, the properties 
of the dust are expected to change at some level, so that Ay /tiqq is 
not only uncertain in the absolute sense but may vary from location 
to location, e.g., due to grain coagulation or the growth of mantles. 
This may produce some of the scatter in our comparison to direct 
extinction measures above. Future observations with ALMA, HST, 
and ground based photometry all offer the potential to improve this 
work substantially. In the meantime, we present a first-order com¬ 
parison using the best publicly available data. 


2.3.5 Comparison Between Herschel and Spitzer Results 

Before longer wavelength Herschel data were available, we mod¬ 
elled IR emission in the LMC using Spitzer 70 and 160 pm maps 
from SAGE survey (Meixner et al. 2006). We used comparisons to 
coarser resolution data at lOO/rm to help us account for the out- 
of-equilibrium emission from very small grains (VSGs) contribut¬ 
ing to the IR emission at 70 pm, finding about 50% of the emis¬ 
sion to represent contamination but otherwise the approach was 
very similar to our main results here. Table 1 compares the me¬ 
dian and standard deviation of Tduat and Ay in the MAGMA field 
between Spitzer and Herschel. On average, Tduat from Spitzer is 
~ 11% lower than our best-fit value from Herschel, while we find 
a ~ 32% higher Ay using Spitzer than we do with Herschel. The 
point-by-point correlation between the maps is good, with a Pear¬ 
son coefficient of p = 0.76 and 0.9 from comparing Tduat and 
Ay between maps derived from the two telescopes. Consequently, 
the qualitative results of this paper would remain unchanged if we 
use dust properties derived from either telescope. For our purposes, 
the main change would be that the slope in the Ico-Ay relation is 
somewhat lower if we use only Spitzer data. 

As a sanity check, we also compare to Skibbaetal. (2012), 
who used the HERITAGE data to derive dust temperature and dust 
mass in the Magellanic Clouds. They fit a modified black body with 
j3 = 1.5 at A 300 pm and a different /3 that is allowed to vary to 
best describe the observation at 350 and 500 pm. We relate our es¬ 
timate of Tieo to their dust mass estimate using their adopted mass 
absorption coefficient, Kieo ~ 13.75 crn^g”^. Using this Kieo, 
Dlmc = 50 kpc, and Dsmc = 60 kpc, our LMC neo map con¬ 
tains ~ 1.0 X 10® Mq, within 10% of the ~ 1.1 x Itf Mq found by 
Skibba et al. (2012) and also similar to the « 1.2 x 10® Mq found 
using Spitzer data by Bernard et al. (2008). In the SMC field, we 
find 0.9 X 10® Mq in our bright lines of sight, again just slightly 
below the ~ 1.1 x 10® Mq from Skibba et al. (2012). Considering 
potential issues with aperture matching and subtleties of fitting, our 
maps appear very consistent with previous works. 
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Figure 2. Dust temperature (left) and Ay (right; from rieo) maps of the LMC (upper panels) and SMC (lower panels). The fields covered by CO data are 
shown in black contours. In addition to the majority of bright CO, the MAGMA field (black contours in the LMC) also shows enhancement of visual extinction 
relative to other regions in the LMC. That is, there is already a “by-eye” IcO'^V relation apparent in the figure. (A colour version of this figure is available 
in the online journal.) 


2.4 Coverage and Av Distribution in the Magellanic Clouds 

As a targeted follow-up survey, MAGMA does not cleanly sample 
the distribution of Av in the LMC. Instead, MAGMA preferen¬ 
tially samples high extinction lines of sight. Likewise, a similar bias 
in Av is expected for the APEX and SEST fields in the SMC. Fig¬ 
ure 2 shows the maps of dust temperature {Tdust) and dust optical 
depth (ri6o) scaled to visual extinction (Av) in the LMC (upper 
panels) and SMC (lower panels) fields. The solid black contours 
show the regions covered by CO maps in each galaxy, the MAGMA 
field for the LMC and the APEX and SEST fields for the SMC. The 


upper-right panel clearly shows that dust shielding (estimated from 
ri6o) in the MAGMA field is enhanced compared to other regions 
in the LMC. This also appears to be the case for the APEX and 
SEST fields in the SMC. Considering the fact that the MAGMA 
field harbors the brightest molecular clouds identified by previous 
CO surveys in the LMC (Fukui et al. 1999, 2008), even this simple 
visual comparison implies a close relation between Av and Ico in 
the Magellanic Clouds. In Figure 2, we also note a quite significant 
variation of dust temperature across the Magellanic Clouds. This 
reflects the variation of interstellar radiation field (ISRF) strength. 
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Figure 3. Top: Histogram of Ay distribution in the LMC field (grey filled 
histogram; it includes all regions in the LMC with the Herschel IR intensi¬ 
ties greater than their 3 -(t uncertainties) and the MAGMA field (red filled 
histogram; this is the region in the LMC field where Ihe MAGMA survey 
mapped in CO). Likewise, the unfilled blue histogram shows distribution 
of Ay in the SMC field, and histogram shaded with blue diagonal lines 
shows that in the CO mapped regions in the SMC (the APEX and SEST 
fields). Low extinction (Ay ^ 1 mag) lines of sight dominate the Magel¬ 
lanic Clouds and also comprise a large portion of the CO mapped area in 
each galaxy. Bottom; Completeness of Ay in the CO mapped regions in the 
Magellanic Clouds (shown as red filled histogram for the MAGMA field, 
and histogram shaded with blue diagonal lines for the APEX and SEST 
fields), calculated by dividing the number of pixels with a given Ay in 
the CO mapped region by the total number of pixels with that Ay in each 
galaxy. The completeness level drops to fa 50 per cent at Ay ss 0.8 mag 
in the MAGMA field and at Ay fa 0.6 mag in the APEX and SEST fields. 
(A colour version of this figure is available in the online journal.) 


which may impact the amount of CO emission in the region via 
photodissociation (e.g. Israel 1997). We further explore this idea in 
Section 3.3, where we divide the MAGMA field into high and low 
Tdust regions. 

In Figure 3, we show histograms of Ay distribution over the 
whole IR-bright area and specifically over the CO mapped regions 
in the LMC and SMC. The bottom panel shows the completeness 
of Ay coverage by the CO map in each galaxy, that is, the fraction 
of LMC and SMC pixels in the specified Ay bin (binsize of 0.1 
dex mag) that lie within each galaxy’s CO map. Therefore, 100% 
in the bottom panel means that all pixels in that Ay bin lie within 
CO map’s field of view. 

Low extinction (Ay ^ 1 mag) lines of sight make up most of 
the area in the Magellanic Clouds, even within the CO mapped re¬ 
gions where we observed enhanced dust shielding relative to other 
regions of the LMC and SMC in Figure 2. On ~ 10 pc scale, the 
distribution of Ay in the MAGMA field is well described by a log¬ 
normal function with mean 0.65 mag and standard deviation 0.30 
dex. The APEX and SEST fields in the SMC shows a narrower 
Ay distribution, which can be fit by a log-normal function with 
mean 0.58 mag and standard deviation 0.17 dex. This is in qual¬ 
itative agreement with the observations of local molecular clouds 


Table 2. IcO'Ay relation fitting parameters in the Pipe nebula and the 
Perseus molecular cloud 


Cloud 

lo 

k 

Akl2 

j^mid 

b 


(K kms-l) 

(mag-l) 

(mag) 

(mag) 


Pipe“ 

32.3 

0.694 

n/a 

4.55 

0.036 

Perseus^ 

42.3 

0.367 

0.580 

n/a 

n/a 


“ /co = 7o((l + . Lombardi et al. 2006 

>> 7co = 7o(l - e-'=(Av-Afci 2 )) . Pineda et af 2008 

on sub-pc scale, where the column density distribution is well fit 
by a log-normal function, often accompanied by a power law tail 
(Kainulainen et al. 2009). At matched spatial resolution, we will 
see that the average extinction through a local Milky Way cloud is 
« 1 — 2 mag (Section 2.5.2), a few times higher than the mean 
Ay we see in the CO surveyed regions in the Magellanic Clouds. 
The Ay associated with a Milky Way cloud would thus represent a 
bright spot, but not a dramatic outlier, in the Magellanic Clouds on 
~ 10 pc scale. 

The completeness of the MAGMA coverage exceeds 50 per 
cent at Ay ~ 0.8 mag. That is, about one half of the pixels in the 
LMC field with Ay ~ 0.8 mag lie within the MAGMA field. In the 
most extreme case where MAGMA recovers all of the CO emission 
from LMC then the bias in our Ico estimate at a given Ay will 
simply be the completeness in that Ay range. So above « 0.8 mag, 
completeness introduces no more than a factor of two uncertainty. 
In reality, MAGMA does not recover all CO emission from the 
LMC and we do not expect the measurement in the MAGMA field 
to be quite so strongly biased. In the SMC, the completeness of the 
APEX and SEST coverage becomes 50 per cent at lower Ay than 
MAGMA, at Ay ~ 0.6, which is expected since the SMC CO map 
covers a large contiguous area in the southwest of the SMC. 


2.5 Milky Way Comparison Data 

We compare the Ico-Ay relation in the Magellanic Clouds to 
the Milky Way using three data sets: (1) analytic approximations 
to the highly resolved (sub-pc) Ico-Ay relation measured in the 
Pipe nebula and Perseus molecular clouds, (2) observations of lo¬ 
cal molecular clouds degraded to ~ 10 pc resolution, and (3) the 
pixel-by-pixel relation for high Galactic latitude lines of sight also 
convolved to ~ 10 pc resolution. 


2.5.1 Highly Resolved Milky Way Clouds 

The proximity of Milky Way molecular clouds allows highly re¬ 
solved comparisons of Ico and Ay, with the limiting reagent 
mostly wide field CO maps. We are aware of two quantita¬ 
tive studies of the dependence of Ico on Ay in nearby clouds: 
Lombardi et al. (2006) considered the Pipe nebula and Pineda et al. 
(2008) studied Perseus. Pineda et al. (2010) carry out a similar 
study of Taurus but do not analyze the relation in exactly the way 
we need for this comparison. These studies find 


Ico 


+ -6) :Pipe 

7o(l — : Perseus 


(3) 


where for Perseus /q is the integrated intensity at saturation, Aki 2 
is the minimum extinction needed to get CO emission, and k is the 
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conversion factor between the amount of extinction and the optical 
depth. In the Pipe, the relation looks similar, but here /o(l — fe) is 
the saturation intensity, and the minimum extinction required for 
CO emission is equals to Av ~ In We list the 

best fit parameters'^ reported in the above studies for these clouds 
in Table 2 and plot the two relationships as a point of comparison 
throughout the paper. 

The qualitative behaviour of the relations observed in the Pipe 
and Perseus highlight some of the key physics governing CO emis¬ 
sion from molecular clouds (e.g., classic PDR models such as 
Maloney & Black 1988; van Dishoeck & Black 1988). First, there 
appears to be a minimum amount of dust extinction required for 
bright CO emission. Below this level, CO abundance is very low 
because of photodissociation, leading to no or negligible CO emis¬ 
sion. Above that threshold there is an approximately linear relation 
between Av and Ico for some range of Av- Then at very high 
Av, CO intensity saturates as the line becomes very optically thick 
across the whole velocity range. In this optically thick regime the 
observed CO intensity becomes a product of excitation temperature 
(Tsx), the beam filling factor, and linewidth. While this observed 
dependence of Ico on Av highlight the importance of dust shield¬ 
ing for CO emission, the differences among the relation for the two 
clouds and even within an individual cloud make it clear that Av 
is not the only factor that determines the amount of CO emission. 
Different geometries and environmental factors (external radiation 
field, density structure, internal heating) will lead to a substantial 
dispersion in CO emission even for the same amount of shielding. 
Indeed, one of the main conclusions of Pineda et al. (2008) was 
that environmental effects can be very strong even within the same 
molecular cloud complex (e.g., see their Figure 6). 

2.5.2 Integrated Measurements for Milky Way Clouds 

To make a more direct comparison of the Milky Way to the Magel¬ 
lanic Clouds, we consider local Galactic molecular clouds from Ta¬ 
ble 1 in Dame et al. (2001). We calculate their average Av and Ico 
values on ~ 10 pc scales to simulate how they might appear in one 
of our Magellanic Cloud maps and report the measurements in Ta¬ 
ble 3. We compile the Av and Ico values of these clouds from the 
Milky Way CO map by Dame et al. (2001) and the E{B — V) map 
published by the Planck collaboration (Planck Collaboration et al. 
2013a), where we take Rv = 3.1 to convert E{B — V) into 
Av- We refer the interested reader to Planck Collaboration et al. 
(2013a) for more detailed information on the Planck E{B — V) 
map and note that Av in local molecular clouds calculated from 
the Planck E{B — V) map may be slightly overestimated ®. 

^ Note that the best fit parameters for the Pipe nebula give negative mini¬ 
mum extinction threshold, which means that Ico positive at Av = 0. 
These authors also fit Ico ~ Av relations for these clouds using a linear 
function (i.e. with the functional form Ico = Avo + rAv, where Avo 
is the minimum extinction threshold for CO emission and r is the linear 
coefficient relating Av and Ico), tmd in this case the minimum extinction 
threshold for CO emission in the Pipe nebula is positive. 

® The Planck E{B — V) map used in this analysis is a conversion of their 
dust optical depth (r) to E{B—V) using the correlation between E(B—V) 
of SDSS quasars and dust optical depth at high Galactic latitude lines of 
sight (i.e., the conversion we discussed in Section 2.3.2, see Figure 22 
in Planck Collaboration et al. 2013a). However, Planck Collaboration et al. 
(2011c) compared column density Af(H) from NIR extinction with dust 
optical depth, finding higher dust emissivity, T/Af(H), in the molecu¬ 
lar phase by a factor ~ 2 than the atomic phase (see also Figure 20 in 


To simulate the ~ 10 pc physical resolution of our Magellanic 
Cloud data, we convolve the Planck Av map (original resolution 
5') and the Dame et al. (2001) CO map (original resolution 7.5') 
with a 3° Gaussian kernel. At 200 pc, a typical distance to a nearby 
molecular cloud, this resolution corresponds to « 10.5 pc. Most 
local clouds subtend several pixels at this resolution (7 pixels for 
the Chamaeleon molecular cloud at the least; at the most 313 pix¬ 
els for the Taurus-Perseus-Auriga complex). We take the median 
Av and Ico across all of these pixels as representative and report 
them in Table 3; we adopt the Av and Ico of each cloud at zb 1- 
cr percentiles as the representative uncertainties when plotting the 
data. 

2-5-3 High Galactic Latitude (\b\ > 5°) Emission 

As a final point of comparison, we use all-sky Milky Way CO and 
Av maps to explore the high Galactic latitude lines of sight. We use 
the same Planck E{B — V) map® described above section to calcu¬ 
late Av- Because the Dame et al. (2001) CO map only covers lim¬ 
ited range of Galactic latitudes (|6| < 30°), we use the the Planck 
TYPE 2 CO map. This is a map of integrated CO line emission that 
has been extracted from the Planck HFI channels using the multi¬ 
channel method by the Planck team (Planck Collaboration et al. 
2013b). This map has an angular resolution of 15', and is known 
to be better suited for intermediate / high Galactic latitude regions 
than the TYPE 1 CO map (Planck Collaboration et al. 2013b). 

Before comparing this map with the Av map, we corrected 
for the contribution of ^®CO to the map by dividing the TYPE 2 
CO map by 1.2. We then degraded the resolution of the CO map 
to match that of the Av map. Using the Av and CO maps de¬ 
scribed above, we construct a pixel-by-pixel Ico-Av relation for 
the high latitude Milky Way sky suitable for comparison to Mag¬ 
ellanic Cloud measurements. We consider all area with |ti| > 5°, 
avoiding the Milky Way midplane in order to remove confusion 
from multiple components along a line of sight. 

3 RESULTS 

3.1 Ico vs. Av in the Magellanic Clouds 

Eigures 4 and 5 show our main observational results, the Ico-Ay 
relationship in the Magellanic Clouds. Eor this analysis, we con¬ 
sider the MAGMA field for the EMC and the APEX field for the 
SMC. Eor the APEX CO data, we assume a ratio of unity to trans¬ 
late CO J = 2 —>• 1 to CO J = 1 —>■ 0 (e.g., see Bolatto et al. 
2003). For clarity, the Ico-Av relationship in the N83 complex is 
not displayed in Figures 4 and 5, but we note that it is similar to the 
Ico-Av relationship in the southwest region of the SMC. 

Figure 4 shows the relationship between Ico and Av pixel- 
by-pixel for each galaxy (the EMC in upper panels and the SMC in 
lower panels). The distribution of individual data points are shown 

Planck Collaboration et al. 2013a). Therefore, we caution that the Planck 
E{B — V) map may systematically overestimate the actual E{B — V) 
(and Av) in local clouds. For example, comparison of E(B — V) from 
NIR extinction and the Planck E{B — V) map in Taurus and Ophiuchus 
molecular clouds (Planck Collaboration et al. 2013a) show that correlations 
between them are pretty strong (their Table 5), but the Planck E(B — V) 
is systematically higher by 25 per cent. 

® Unlike the case for the local clouds, high Galactic latitude lines of sight 
are diffuse and the application of E{B — V)/r conversion is not expected 
to overestimate the actual E{B — V) value in this case. 
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Figure 4. Pixel-by-pixel comparison of Ay (x-axis) and Iqq (y-axis) in the LMC (top panels) and SMC (bottom panels), shown in semi-log scale in the left 
panels and in log scale in the right panels. In the left panels, the distribution of data points in the LMC and SMC are shown as the contours, while the grey 
points show individual lines of sight where the density falls below 10 points per bin. The lowest contour means there are at least 10 individual data points 
in the bin, increasing by factors of 2 with each step for the LMC and by factors of 3 for the SMC. The red circles with en'or bars indicate the median CO 
intensity and its scatter measured in bins of log Ay, our LMC and SMC Ico-Ay relations. For comparison, we also plot the I^Q-Ay relations measured 
at much higher (~ sub-pc) resolution for two Galactic clouds, Perseus (dash-dotted, Pineda et al. 2008) and the Pipe (dashed : Lombardi et al. 2006). Blue 
triangles represent a more well-matched Milky Way comparison, Galactic moleculai' cloud data compiled from Planck (Planck Collaboration et al. 2013a) and 
Dame et al. (2001), degraded to physical resolution similar to what we achieve in the Magellanic Clouds (3°, which con'esponds to 10.5 pc at the distance 
of 200 pc). In the semi-log plots one can see how median CO intensity approach zero as Ay decreases in the Magellanic Clouds, while in the log-log plots we 
rescale the y-axis to better compare I^Q-Ay relations in the Magellanic Clouds with Galactic moleculai* clouds. The horizontal eiTor bar in the bottom-right 
comer of each panel shows a typical uncertainty in Ay derived from Monte-Carlo analysis (Section 2.3.3). (A colour version of this figure is available in the 
online journal.) 
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Table 3. Median Ay and Iqq of Galactic Molecular Clouds at 3° resolution 


Cloud 

(degrees) 

ba 

(degrees) 

Area“ 

(deg2) 

Pixels 

< Ay >*’ 
(mag) 

< Ico 
(K km s“^) 

< Ico/Ay 
(K kms“^ mag~^) 

Aquila 

20 

8 

227 

85 

4.0 

7.5 

2.0 

Camelopardalis 

148 

20 

159 

61 

0.54 

0.32 

0.55 

Chamaeleon 

300 

-16 

27 

7 

1.4 

2.0 

1.6 

Gum Nebula 

266 

-10 

97 

37 

1.5 

0.46 

0.31 

Ophiuchus 

355 

17 

422 

151 

1.9 

2.6 

1.3 

Orion 

212 

-9 

443 

163 

2.1 

1.6 

0.65 

Polaris Flare 

123 

24 

134 

55 

0.81 

2.5 

2.7 

Taurus 

170 

-15 

883 

313 

1.9 

3.5 

1.8 


“ Compiled from Table 1 in Dame et al. (2001). 

^ Computed from Planck E{B — V) map, assuming Ry = 3.1. Note that this Ay may be systematically overestimated (see text). 

Computed from Dame et al. (2001) CO map. 
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Figure 5. Average Iqq-Av relation in the LMC (grey downward triangles) 
and SMC (red circles) from Monte-Carlo simulation (Section 2.3.3). Other 
lines and symbols are the same as Figure 4. Here, the vertical en'or bars as¬ 
sociated with the average relations do represent uncertainty in the average 
relations. Note that high Ay points in the average relations arise from very 
few lines of sight moving across Ay bins as we simulate Ay maps and 
thus mimic saturation effect observed in highly resolved Galactic molec¬ 
ular clouds (see Section 2.5.1). The agreement between the LMC, SMC, 
and integrated measurements for the Milky Way is fairly good at Ay ^ 4, 
suggesting that dust shielding is the primary factor that determines the dis¬ 
tribution of bright CO emission from a molecular cloud. (A colour version 
of this figure is available in the online journal.) 


as the contours and grey points in the left panels (semi-log scale). 
The median and 1-a scatter of the data are shown as the red circles 
(binned by log Ay) in both the left (semi-log scale) and right (log 
scale) panels. In the left panel of Figure 4, we find that the 7co- 
Ay relation for individual lines of sight shows large scatter, but 
the median trend (red circles) confirms the impression from simple 
visual comparison in Figure 2. That is, the CO intensity appears to 
increase as a function of Ay in the Magellanic Clouds. 

In Figure 5, we plot this average relation for both galaxies. 


Here the mean and uncertainty derive from Monte-Carlo simula¬ 
tion (Section 2.3.3), so that the error bars reflect the uncertainty 
in the mean across 1,000 realizations for the LMC and an equiva¬ 
lent uncertainty for the SMC. These average relations in Figure 5 
are thus our best estimate of how CO intensity (Ico) depends, on 
average, on dust shielding (Ay) at cloud-scales in the LMC and 
SMC. The estimate takes into account the systematic uncertainties 
in estimating Ay from the IR emission. Interpreting the error bars 
requires some care; many of the factors that we simulate create 
correlated errors across a whole realization. Therefore one should 
largely view the error bars as the space within which the mean re¬ 
lation could move. 

Care must also be taken when interpreting the high Ay por¬ 
tion of the average relations. Recall that low Ay (Ay ^ 1) dom¬ 
inate the Magellanic Clouds, and there are very few data points 
in the high Ay bins in the LMC and SMC (e.g. see Figure 3). 
For example, there are only 15 lines of sight that has Ay greater 
than 10 in the LMC, and no line of sight has Ay greater than 4 
in the SMC. The greater number of points at lower Ay will lead 
those points to preferentially scatter to high Ay and contaminate 
the measurements at high Ay in the Monte Carlo calculation. This 
will artificially lower the mean Ico at high Ay in the Monte Carlo 
simulation. 

With this in mind, we note that except at Ay 4 in Figure 5, 
the LMC and SMC in Figures 4 and 5 overlap one another, showing 
similar CO intensity at a given Ay. This agreement between the 
LMC and SMC despite their factor of « 2 difference in metallicity 
is consistent with the theoretical picture that CO intensity depends 
on dust shielding (Ay) in an approximately universal way. In the 
next Section, we will compare the Magellanic Clouds to the Milky 
Way in several ways to see that this holds true from Z « 1/5 Zo 
up through solar metallicity. 

We noted earlier that the measured Jco scatters considerably 
at a given Ay in Figure 4. Much of this scatter reflects the noise 
in the CO map, and this scatter clearly dominates the distribution 
of individual data below Ay « 1 mag. Still, we verified that a real 
CO signal emerges as we stack the CO spectra in the LMC in bins 
of Ay (see Figure Bl). At higher extinctions. Ay ^ 1 mag, we 
also see substantial additional scatter about the median CO inten¬ 
sity in fixed Ay bins. We interpret this as real astronomical signal, 
indicating that line-of-sight Ay on 10 pc scale is not a perfect pre¬ 
dictor of Ico even in the absence of noise. This is reasonable and 
expected as many physical effects beyond shielding may influence 
CO emission, for example, variations in cloud structure, geometry. 
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chemistry, and the interstellar radiation held. The important caveat 
to bear in mind is that while Ico does depend on Ay, the relation¬ 
ship emerges only after substantial averaging because individual 
lines of sight have large scatter in Ico-Ay parameter space. 


3.2 Comparison to the Milky Way 

We are primarily interested in testing the hypothesis that CO emis¬ 
sion depends on dust shielding (Ay) in the same manner across 
environment that differ in metallicity. In the previous Section, we 
saw that CO intensity at a given Ay is indeed similar in the LMC 
and SMC. In this Section, we compare our Magellanic Cloud re¬ 
sults to the Milky Way in three different ways. 

First, we compare the Magellanic Clouds to the highly re¬ 
solved (sub-pc resolution) Ico-Ay relation measured in the Pipe 
nebula and Perseus molecular clouds. These appear as the dashed 
line (Pipe) and the dash-dotted line (Perseus) in Figures 4 and 5. 
Overall, the Ay-Ico relation in the Magellanic Clouds resemble 
that in the Pipe nebula (note that the Pipe curve is not a fit to our 
data or even normalized to match our data), while the Magellanic 
Clouds exhibit 2-3 times fainter CO emission than the Perseus 
molecular cloud at high Ay. Most importantly, neither the LMC 
nor the SMC in Figure 4 exhibit the qualitative features observed 
in the Pipe and Perseus (Section 2.5.1), at least not prominently. 
We observe no clear saturation in Ico and the threshold behavior, 
if present, appears far weaker than in the Perseus molecular cloud. 
The average Ico-Ay relations in Figure 5 seem to exhibit a satu¬ 
ration of CO at high Ay (Ay ^ 10 mag for the LMC, Ay ^ 4 
for the SMC), but as noted earlier, the apparent trends observed at 
high Ay in the average Ico-Ay relations are vulnerable to arti¬ 
facts arising from Monte-Carlo simulation. 

The lack of evidences for the saturation of CO and threshold 
behaviour in the Magellanic Clouds almost certainly stem from the 
dramatic difference in resolution between our measurement in the 
Magellanic Clouds and those used to construct the Pipe and Perseus 
curves. Our measurements combine large parts of a cloud (~ 10 
pc) into a single data point, so that each line of sight is an average 
of Ay and Ico in a 10 pc beam. On the other hand, the Pipe and 
Perseus relations are measured at much higher resolution (sub-pc). 
These clouds would not exhibit Ay ~ 5-10 mag at 10 pc resolu¬ 
tion. For a fairer comparison to these highly resolved curves, one 
would need very high resolution Ay and Ico data in the Magel¬ 
lanic Clouds. 

An alternative approach is to measure Ay and Ico averaged 
over a 10 pc beam for the local clouds. These appear as blue trian¬ 
gles in Figures 4 and 5. These triangles represent our best estimate 
of how local clouds would look like at the distance of the Mag¬ 
ellanic Clouds. This makes them analogous to the individual grey 
points in Figure 4. At matched ~ 10 pc resolution, most of the 
Milky Way clouds have rather low Ay, ranging from ~ 0.5 mag to 
~ 4.0 mag. The agreement between the Milky Way and the Magel¬ 
lanic Clouds is much better in the case of matched resolution than 
for the highly resolved Milky Way relations. A close inspection of 
Figure 5 suggests that some of the local clouds have rather fainter 
CO emission at a given Ay than the average lines of sight in the 
Magellanic Clouds, but this could be partly due to overestimated 
Ay in the local clouds in the Planck maps (Section 2.5.2). Given 
that the scatter in Ico is very large at a given Ay both for the 
Magellanic Clouds and the Milky Way clouds, we interpret the in¬ 
tegrated measurements of Ico-Ay in local clouds to substantially 
agree with our measurements in the Magellanic Clouds. This rein- 
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Figure 6. Integrated CO intensity (y-axis) as a function of Ay (at-axis) at 
high Galactic latitude (|t)| > 5°) lines of sight in the Milky Way, compiled 
from Planck data (Planck Collaboration et al. 2013a) degraded to ~ 3° res¬ 
olution. The contours and grey points show the distribution of data points in 
the Ico-Ay space, following the same contour definition for the LMC in 
the left panel of Figure 4. The median CO intensity at a given Ay is shown 
as the black triangles, and the average LMC Ico-Ay relation (Figure 5) 
is shown as the red circles, shifted slightly toward right in x-axis for eas¬ 
ier comparison. (A colour version of this figure is available in the online 
journal.) 

forces our results from the internal comparison of the two Magel¬ 
lanic Clouds. 

As a final comparison, we plot the Ico-Ay relation for high 
latitude (|&| > 5°) Milky Way lines of sight in Figure 6. The 
high latitude Milky Way lines of sight (black triangles) exhibit a 
similar shape to the Magellanic Cloud relations but with a ten¬ 
dency towards lower Ico at a given Ay. That is, the Magel¬ 
lanic Clouds appear to be brighter than the high latitude Milky 
Way in CO on average for a given Ay. Quantitatively, the me¬ 
dian Ico!Ay for Milky Way high latitude lines of sight is « 
0.5 K km s“^ (mag)“^, which is about 0.25 times the value for 
the LMC (« 1.9 K km s"^ (mag)”^). 

Our best explanation for this difference is that the Milky 
Way values are likely to be biased low by dust associated with a 
long path length through the Milky Way Hi disk. For example, at 
b ~ 10° a 200 pc thick Hi disk will yield an integrated path length 
of roughly a kpc. This path length will preferentially sample atomic 
gas, which has a higher scale height than molecular gas and lit¬ 
tle associated CO emission. More, the dust along that line of sight 
through and extended disk will contribute little to shielding distant 
CO from dissociating radiation. Such effects will undermine any 
mapping between line of sight extinction and the local dust shield¬ 
ing that should affect CO emission. In Appendix A2, we show that 
a simple correction for dust associated with an extended Hi disk 
leads to a median Ico at a given Ay for the high latitude Milky 
Way that more closely resembles what we find in the Magellanic 
Clouds. 

The Magellanic Cloud data may also be biased high by our 
focus on the CO mapped regions. As shown in the bottom panel 
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of Figure 3, the completeness of Ay in the MAGMA field drops 
to ~ 50 per cent at Ay ~ 0.8 mag, and at Ay ^ 0.6 mag in 
the APEX and SEST fields. This means that at Ay ~ 0.8 mag, 
IcolAy in the EMC has a lower limit of about half the current 
value if we assume that the other half lines of sight not mapped 
in CO do not have associated CO emission at all, and a similar 
logic can be applied to the SMC as well. If this is the case, Ico at a 
given Ay would be closer between the Magellanic Clouds and high 
latitude Milky Way. Even if this bias drives the results. Figure 6 
suggests the somewhat surprising result that the active parts of the 
Magellanic Clouds are better at emitting CO than the high latitude 
Milky Way. 

Overall, the sense of the comparison made in Section 3.2 is 
this: the Magellanic Clouds, at least the bright regions covered by 
CO surveys, show Ico at a given Ay comparable to or somewhat 
below than those found at highly resolved (sub-pc) local Galactic 
clouds. The Magellanic Clouds data do not show clear evidence for 
saturation of CO line and Ay threshold observed in high resolution 
Milky Way data, likely due to low resolution. After accounting for 
resolution differences, we find that a sample of Milky Way clouds 
at 10 pc resolution largely overlaps the average Ico-Ay relations 
in the Magellanic Clouds. Taking a broader view and considering 
high latitude emission from the Milky Way, we find a qualitatively 
similar Ico-Ay relation to the Magellanic Clouds but note impor¬ 
tant quantitative differences with the sense that at a fixed Ay, gas 
in the high latitude Milky Way emits somewhat less CO than the 
regions covered by the CO maps in the Magellanic Clouds. 

3.3 Influence of the Interstellar Radiation Field 

Figure 4 showed large scatter in Ico at fixed Ay, implying addi¬ 
tional physics beyond the abundance of dust at play. One simple 
and often-discussed “second parameter” is the interstellar radiation 
field (ISRF Israel 1997; Pineda et al. 2009). We have estimated this 
quantity as part of our dust modelling: the dust temperature, Tdust, 
depends on the strength of ISRF such that Xisrf ~ 
radiation field drives the photodissociation of CO molecules, so that 
we may reasonably expect Tdust to affect Ico within a bin of fixed 
Ay. 

Figure 7 shows a simple test for the effect of Tdust as a second 
parameter (as an extension to this test, we also discuss the effect of 
strong radiation field on the Ico-Ay relation in Appendix A3, by 
focusing on 30 Doradus complex in the EMC). For this experiment, 
we divide our EMC data into two bins: a “cool” region with lower- 
than-median Tdust, i.e., Tdust < 24.1 K, and a "hot” region with 
higher-than-median Tdust, Tdust > 24.1 K (See Figure 2 for the 
dust temperature map). The median dust temperatures of these two 
regions are ~ 26.0 K (“hot”) and ~ 22.7 K (“cool”). With our fidu¬ 
cial /3 = 1.5, this temperature difference corresponds to difference 
in ISRF strength of a factor of ~ 2.1. Note that the ISRF in the 
MAGMA field is already on average greater by a factor of ~ 5.8 
than the Solar Neighborhood, where Tdust ~ 17.5 K. We calculate 
the average Ico-Ay relations in the “hot” and “cool” regions and 
estimate their uncertainties using Monte-Carlo simulation. 

Interestingly, Figure 7 shows that Ico at fixed Ay is higher 
in the cool region (blue downward triangles) than the hot region 
(red upward triangles) in the MAGMA field. This observed trend 
supports the idea that an intense ISRF suppresses CO emission, 
since several effects might lead to the opposite direction of what 
we see in Figure 7. First, the uncertainties on neo and Tdust are 
correlated in the sense that underestimating Tdust would lead us to 
overestimate rieo. This has the opposite sense of what we see here 
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Figure 7. Temperature dependence of EMC Ico-Ay relation. We divide 
the MAGMA field into hot (Tdust > 24.1 K; the red upward triangles with 
error bars, shifted slightly toward right in x-axis for easier comparison) and 
cool (Tdust < 24.1 K; the blue downward triangles with error bars) regions 
based on the median dust temperature ('--' 24.1 K). Interestingly, the region 
with hotter dust temperature (and more intense ISRF) shows fainter CO 
emission measured at a fixed dust shielding (Ay) compared to the region 
with colder dust temperature. This is the opposite of what one would expect 
for the region with hotter dust temperature having stronger CO emission 
due to a higher kinetic temperature, and we interpret this as the evidence 
for an intense ISRF to suppress CO emission at a fixed dust shielding in 
the hot region more than the cool region. (A colour version of this figure is 
available in the online journal.) 

and so seems unlikely to drive the separation in Figure 7. Second, 
if a higher ISRF leads to higher excitation for the CO, we might 
expect a higher brightness temperature in regions with high Tdust- 
Again, this has the opposite sense of the separation we observe. A 
correlation between Tdust and (3 might produce some of the sig¬ 
nal we see, but such correlations are notoriously difficult to verify 
(Shetty et al. 2011). The simplest explanation of Figure 7 is that an 
intense ISRF tends to suppress CO emission at a fixed dust abun¬ 
dance, presumably due to enhanced dissociation of CO. Enhanced 
dissociation of CO under a strong ISRF then would lead to reduced 
area of CO photosphere, likely causing more beam dilution and 
decreased integrated CO intensity. This agrees with the results of 
Israel (1997), who use dust emission as a tracer of H 2 column den¬ 
sity to calculate Xco, and find that Xco increases as a function 
of the strength of interstellar radiation field. On the other hand, this 
disagrees with studies focused on GMC properties measured from 
CO emission (Hughes et al. 2010, Pineda et al. 2009), which find 
no apparent dependence of the properties of CO clumps or Xco 
on the strength of radiation field. 

3.4 Median Ico Mv across different environments 

As a summary of our main results, we compile the median 
Ico /Ay values across different environments at matched reso¬ 
lution in Table 4. In the MAGMA field, the median Ico /Ay 
is ~ 1.9 K km s“^ (mag)“^. We find slightly lower values in 
























Table 4. Median Iqq I-^V values across environments 


Galaxy 

Region 

< Ico/^V >“ 

(K kms“^ mag~^) 

LMC 

MAGMA field 

1.9 


Hot 

1.5 


Cool 

2.2 

Milky Way 

High latitude, without Hi correction 

0.5 


High latitude, with Hi connection 

0.8 


Local clouds from Planck 

1.3 


Standard Xqq 

4.7 

SMC 

APEX field (CO J = 2 1) 

1.1 


SEST field (CO J = 2 ^ 1) 

1.5 


SEST field (CO J = 1 0) 

1.4 


“ We caution that these values are very sensitive to the choice of pixels for 
computing the values, baseline correction for CO emission in the case of 
the LMC, and uncertainties associated with deriving Ay from IR emission. 

the APEX field and the SEST field in the SMC, with the me¬ 
dian Ico/Ay varies from 1.1 to 1.5 K km (mag)“^ depend¬ 
ing on which transition of CO line to use in the calculation. The 
median Ico/Ay value for local clouds at matched resolution is 
« 1.3 K km s“^ (mag)“^. The lower Ico/Ay in the integrated 
measurements of local clouds partly arises from overestimation of 
Ay in these clouds in the Planck data (see Section 2.5.2). Also, 
we note that calculating a median value is sensitive to the area in 
which the calculation is made and that cloud boundaries are not 
known with precision. If our adopted cloud areas (Table 3) are too 
large then they may bias the local cloud measurements somewhat 
low. As discussed in Section 3.2, the median Ico/Ay in Milky 
Way high latitude lines of sight is significantly lower than the LMC 
and SMC, but a correction for contamination by an extended Hi 
disk (see Appendix A2 for details) bring these values into closer 
agreement. 

For comparison, we also recast Galactic CO-to-H 2 conversion 
factor (Xco = 2 x 10^'’ cm“^ (Kkm s“^)“^, Bolatto et al. 2013) 
into these units. One can calculate the corresponding Ico/Ay for 
a standard Xco, ~ 4.7 K km s“^ (mag)“^, by adopting a Galac¬ 
tic dust-to-gas ratio (W(H) = 5.8 x 10^^ E(B — V), Bohlin et al. 
1978) and Ry — 3.1. Unlike the rest of the median Ico/Ay val¬ 
ues, the recast of Xco only considers molecular gas in the denom¬ 
inator and therefore represents an upper limit to the number that 
would be measured for a real could. 


4 DISCUSSION 

Our comparison of the Magellanic Clouds to the Milky Way ar¬ 
gues that on 10 pc scales CO intensity tracks dust column den¬ 
sity, expressed as Ay. A logical corollary, though one we cannot 
prove with existing data, is that perhaps Ico depends on Ay within 
clouds in an approximately universal way, at least when averaged 
over a sizable population of clouds. If this is the case, then the in¬ 
tegrated CO emission from clouds, and by extension the CO-to-H 2 
conversion factor, may be thought of as a problem with several sep¬ 
arable parts. First, the distribution of Ay within a cloud will depend 
on the distribution of gas surface densities combined with the dust 
to gas ratio. Second, the CO emission from the cloud will depend 
on the Ay distribution. In the case where we are interested in the 
CO-to-H 2 conversion factor, the dust-to-gas ratio and surface den- 
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Figure 8. Illustration of our PDF-based calculation of Xco • The grey filled 
histogram shows the same gas column density PDF at a range of dust-to- 
gas ratios. The hatched histograms indicate the implied distributions of CO 
intensity for each of these Ay distributions given the Pipe (red diagonal 
pattern) and Perseus (green horizontal pattern) Iqo-Ay relations. The light 
red cross hatched region at Ay ^ 0.3 shows where the transition between 
H2 and Hi will be important. The qualitative difference between Perseus 
and the Pipe is driven by the threshold in the Perseus I^Q-Ay relation, 
which we illustrate with a green dashed line. (A colour version of this figure 
is available in the online journal.) 


sity PDF will also determine what part of the cloud is Hi and what 
part of the cloud is H 2 . In this simplified view, the CO emission 
emerging from a whole cloud can be expressed as 

Ico = f [DGR X PDF {Y.h)) , ( 4 ) 

where Tih is the column density of hydrogen gas; PDF {Tih) 
is the distribution function for gas column densities in the cloud; 
DGR is the dust-to-gas ratio, which translated a value of Eh into 
Ay\ and function that translates a value of Ay into 

a CO intensity. Here the integral is over all values of Eh and PDF 
should formally be the total column density per differential column 
density bin. 

If one is interested in the CO-to-H 2 conversion factor rather 
than the total CO intensity, then the 

^ ^ /eh {DGR,^h) 

Xh2 /H2 {T.h,DGR) pdf (Eh) ’ ^ ^ 

Here the top integral is as above. The bottom integral sums the total 
H 2 mass in the cloud and includes the function (Eh, DGR), 
which indicates whether at the specified DGR gas of column den¬ 
sity Eh is atomic or molecular. 







































14 Cheoljong Lee et al. 


This very simplified view breaks the CO-to-H 2 conversion 
factor into four parts: the distribution of gas column densities in 
a cloud, the dust-to-gas ratio, the relationship between CO and dust 
shielding, and the dependence of the H 2 /H 1 transition on Av. The 
benefit of stating the problem in this simplified way is that each of 
these topics has been studied independently (see the Introduction). 
The last few years have seen substantial work characterizing the 
column density distribution of local clouds. Theoretical work has 
also established a lognormal distribution as the baseline expectation 
for a turbulent cloud. Substantial recent theoretical work has also 
gone in to understanding the H 1 -H 2 transition. Finally, both the¬ 
oretical and observational efforts (including Section 3) have been 
made to understand the dependence of Ico and Av- That is, this 
approach breaks the topic of CO emission from clouds in galaxies 
into a separable problem whose parts may be more tractable than 
the topic considered as a whole. 


4.1 Implications for the Metallicity Dependence of Xco 

This sketch can be used to make an empirically-driven prediction 
for the metallicity dependence of the CO-to-H 2 conversion factor, 
Xco- If we consider the column density distribution of gas to be 
universal, then the distribution of Av is the dust-to-gas ratio times 
this column density distribution. This Av distribution predicts the 
emergent CO emission from the cloud. The conversion factor de¬ 
pends on the ratio of this CO emission to the sum of the distribution 
of gas over the part of the cloud that is molecular. Then by vary¬ 
ing the dust-to-gas ratio and repeating the calculation we can derive 
how Xco changes as a function of DGR in this simple cartoon. 

Schematically, Figure 8 shows this approach. In this cartoon, 
gas obeys a universal PDF (here a lognormal), which is scaled by 
a dust to gas ratio to yield a PDF of Av values for each cloud. 
That PDF appears as a grey normalized histogram. Applying a Av- 
based prediction for Ico, one arrives at a prediction for the CO 
intensity. That appears in red and green here for two such functions: 
the Pipe and the Perseus relation. This emission would be summed 
to get the integrated emission from the cloud. Finally, some part of 
the cloud is atomic (shown by the light red region) and that is not 
book-kept in the sum of the molecular mass. 

To carry out the calculation quantitatively, we use the parame¬ 
terization of Av PDFs for local clouds at solar metallicity listed in 
Table 1 of Kainulainen et al. (2009). In this analysis, we only con¬ 
sider the log-normal part of the Av PDFs (i.e., ignoring the power- 
law tail). We also ignore Av ^ 1 mag in the solar-metallicity 
clouds because the uncertainty in the extinction mapping used to 
derive the cloud PDFs becomes substantial below Av ~ 1 mag 
(though note that the cartoon in Figure 8 does show lower Av in 
the top panel; this is just illustrative). These lognormal distributions 
clipped at Av > 1 mag are our baseline gas distributions. That is, 
we consider the gas column density PDF to be: 


PDF{EhAZq)) 


PDF{AvAZq)) 

DGR{Zq) 


( 6 ) 


where the subscript i refers to one of the Kainulainen et al. (2009) 
Av PDFs and DGR{Zq) is the solar metallicity dust-to-gas ra¬ 
tio. Because we make a relative calculation of Xco. the numerical 
value of DGR{Zq) will cancel out of our results. 

Without specifying what, precisely, DGR is for the Milky 
Way, we can scale these Milky Way PDFs to those we would expect 
for otherwise identical clouds at some fraction of solar metallicity 
by simply dividing values of Av by the relative dust-to-gas ratio. 


That is, by dividing all Av values by 2, we can shift the PDF to 
represent an otherwise identical cloud at half solar DGR- That is, 
we hold PDF (Sir,;) fixed across metallicity and derive Av.i (Z) 
for some sub-solar metallicity, Z, via: 

PDF{AvAZ)) = ■§§i^PDFiAvAZe)) « -^^PDF (AvAZq)) ■ 
® ® (7) 

In the last step, we take Z and DGR to vary linearly with 
one another, but a more complicated dependence (e.g., see 
Remy-Ruyer et al. 2014) could easily be introduced into the for¬ 
mulae. Also note that at Z < Zq, we will include Av < 1 mag in 
our calculation; the uncertainty surrounding low Av is in the de¬ 
termination of the Milky Way cloud PDFs, not in their inclusion in 
the calculation. 

Next, we translate each Av PDF into a predicted CO inten¬ 
sity. To do, we input the Av distribution to one of the Ico-Av 
relationships discussed in the first part of the paper. In this calcu¬ 
lation, we consider four relations: the Pipe, Perseus, the LMC, and 
the high latitude Milky Way. For the former two we use the best-fit 
functions in Table 2. For the latter two we fit polynomial functions 
to the average Ico-Xv relations for the LMC and Milky Way high 
latitude lines of sight (red circles and black triangles with error bars 
in Figure 6, respectively). 

Thus we have a set of realistic PDFs, denoted by subscript 
i, and four potential Ico-Av relations, which we denote with the 
subscript j- We calculate a plausible set of CO intensities emergent 
from each cloud plus relation (i, j) for each of a range of metallic- 
ities, Z: 


PDF (Ico,i.AZ)) = ff°’'^^{AvAZ)) , ( 8 ) 

where PDF {Ico,i,j{Z)) is a PDF of CO intensity at metallicity 
Z for Av PDF i and assumed Ico-Av relation j. 

To estimate Xco, we need to compare the emergent CO in¬ 
tensity to the amount of molecular hydrogen gas, 'En 2 ,i{Z)- This 
requires one additional step, which is to differentiate between H 2 
and Hi in the PDF. Krumholz et al. (2009) and McKee & Krumholz 
(2010) argue that the layer of Hi in a molecular cloud can be 
described as a shielding layer of nearly fixed extinction and 
observations of Perseus support this (Lee et al. 2012) (see also 
Wolfire et al. 2010; Sternberg et al. 2014). In detail, however, the 
depth of this layer may vary with metallicity, the radiation field, and 
volume density (e.g., see Wolfire et al. 2010). Here we will adopt 
a simple approach and adopt a constant Av^ ~ 0.2 mag for each 
cloud. We take all gas with Av below this value to arise purely 
from atomic gas and so do not book keep it when summing the 
PDF to obtain an H 2 gas mass. A more realistic treatment of this 
H 2 -H 1 transition represents a logical extension of this calculation. 
Here we will only note when this becomes a dominant considera¬ 
tion. 

In this approach 


PDT'(Eh2) 


PDF{^h) ifAv>A^^ 
0 ifAv<Af^ 


(9) 


where A™ is the extinction depth of the Hi shielding layer around 
the cloud. 

Combining these equations, we generate Ico and H 2 distri¬ 
butions as a function of metallicity, cloud structure, and Ico-Av 
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Figure 9. Xqq as a function of dust-to-gas ratio (and thus metallicity), estimated from the Ay probability distribution function (PDF) observed for Milky 
Way clouds (Kainulainen et al. 2009) and various Iqo-Ay relations (Pipe {left top) : Lombardi et at 2006, Perseus {right top) : Pineda et al. 2008, LMC {left 
bottom), MW high latitude {right bottom) : polynomial fits to their average Iqq-Ay relations). We compute Xco a function of dust-to-gas ratio following 
the illustration in Figure 8 (see text for details). Each panel shows results for a different functional form of 7co(^V")’ while the different red curves in each 
panel represent the different PDFs from Kainulainen et al. (2009). The dashed vertical line denotes a dust-to-gas ratio below which the assumed Ay threshold 
for the Hi envelope {Ay^) and the power law tail in the Ay PDF become critical in the calculation of CO emission. Thick light grey lines indicate Xcq 
as a function of dust-to-gas ratio with different power law slopes (Xco cx; DGR““ with a = 1, 2 and 3). Just below solar metallicity, Xqq is inversely 
con'elated with metallicity with tractable scatter (q is between 1 and 2). However, by 1/3 solar metallicity the input PDF, lQo{Ay) and assumed Hi layer 
dominate the calculation that Xco scatters by two orders of magnitude. Fundamentally this reflects the instability of using only a very small fraction of the 
cloud area to trace its total mass and the calculation undermines the utility of CO to study analogs of local clouds at even moderately sub solar metallicities. 
(A colour version of this figure is available in the online journal.) 























































































































































































































































































16 Cheoljong Lee et al. 


mapping. Then, the CO-to-H 2 conversion factor is simply the total 
amount of H 2 column density divided by the total amount of Ico'- 

Xco,i,j{Z) ^ JPDFiIco,iAZQ)) 

Xco,i,j{Z@) j PDF {Ico,i,j{Z)) J PDF {En.AZe)) 

( 10 ) 

By calculating Xco relative to the solar metallicity value, the solar 
metallicity DGR drops out. 

Figure 9 shows the result of the procedure illustrated in Figure 

8, Xco as a function of metallicity (dust-to-gas ratio). In Figure 

9, each panel shows the result for the different adopted Ico-Av 
relation. In each panel, different curves indicate different Av PDFs 
drawn from Kainulainen et al. (2009). Gray lines indicate power 
law dependences of Xco on metallicity for comparison. 

This figure illustrates a few points. First, for most adopted 
PDFs and scalings, between about solar metallicity and ~ 1/3 Zq 
we expect Xco oc Z ^-Z that is, we expect a moderately 
non-linear scaling in this regime. This is consistent with a num¬ 
ber of theoretical and empirical results summarized in Bolatto et al. 
(2013), including results from dust (Israel 1997; Leroy et al. 
2011) and star formation scaling arguments (Genzel et al. 2012; 
Schruba et al. 2012; Blanc et al. 2013). It is steeper than most 
virial mass-based results (Wilson 1995; Rosolowsky et al. 2003; 
Leroy et al. 2006; Bolatto et al. 2008). The constraint here can be 
phrased as follows: if local cloud PDFs are rescaled lower dust- 
to-gas ratios with no other change in the cloud physics, we might 
expect Xco on the scale of whole clouds to scale as Z~^ to Z~ . 

Second, the calculation becomes very sensitive to the H 2 -H 1 
prescription, the shape of the PDF, and the adopted Ico-Av rela¬ 
tion at low metallicity. Even as high as 1/3 Zq these factors cre¬ 
ate a substantial (factor of 2) spread among our results. Below this 
value they dominate the results. At some level the spread in our 
estimates corresponds to a spread in nature, so that this highlights 
intrinsic scatter or uncertainty in the use of CO to trace H 2 at low 
metallicities. In this regime, large parts of a cloud may be Hi, any 
threshold for CO emission will become incredibly important, and 
any power law tail or extension to high Ay will be preferentially 
very good at emitting CO. This simple PDF-based approach argues 
for intrinsic uncertainties of order a factor of a few when using CO 
to trace H 2 below about 1/3 solar metallicity. 

Finally, though not our focus, the extension of these trends 
to super-solar metallicity suggests that a factor of ~ 2 change 
in Xco could be achieved by bringing the Solar Neighborhood 
clouds to even higher metallicity (and so stronger shielding). At 
the same time one might expect a number of other changes in 
the ISM, such as the emergence of a widespread diffuse molecu¬ 
lar phase. But put simply, if all of the molecular gas in the Milky 
Way were better shielded, as one might expect for identical clouds 
dropped into a more dust-rich system, the conversion factor might 
be expected to be a factor of ~ 2 lower (Planck Collaboration et al. 
2011b). This is not clearly observed (e.g., Donovan Meyer et al. 
2013; Sandstrom et al. 2013), but is also not clearly ruled out by ob¬ 
servations given that the link between Milky Way and extragalactic 
observations remains uncertain at the 10s of percent level. 


4.2 Physics, Key Unknowns, and Complicating Factors in the 
Ico-Av Relation 

In calculating Xco, the Perseus Ico-Av stands out because it in¬ 
cludes a hard Ay threshold for CO emission. This is particularly 
stark in Figure 8 as below about 1/3 solar metallicity almost none 


of the PDF exceeds the Perseus threshold, suggesting an almost 
totally CO-dark cloud. Such a threshold is not visible in our Mag¬ 
ellanic Cloud measurements, nor is the saturation in CO intensities 
found at high Ay in both the Pipe and Perseus. This is more likely 
a reflection of our coarse resolution than the absence of these phys¬ 
ical features in the Magellanic Clouds. 

In the Appendix, we demonstrate the presence of substantial 
beam dilution in stacked LMC spectra by comparing line widths — 
which stay about constant — and peak temperatures — which drop 
to unphysical low levels at low Ay. This means that we read our 
observations results as consistent with a universal sub-resolution 
Ico-Ay relation, but not as proof of such a relation. Future mea¬ 
surements comparing dust column density to CO emission across 
diverse environments will be needed to establish this relation and 
its variation at high (~ 0.1 pc, matched to molecular cloud sub¬ 
structure) resolution. Doing so, key questions will be: 

(i) What is the form of the Ico-Ay relation at low extinction? Is 
there a threshold or steepening of the relation at « 0.5 mag? Does 
this also appear in highly resolved maps of the Magellanic Clouds? 

(ii) How does the Ico-Ay change with ambient radiation field? 
In a first study of this sort, Indebetouw et al. (2013) found evidence 
for suppressed CO emission at low Ay in regions illuminated by 
the strong radiation field of 30 Doradus. 

(iii) Is the saturated regime important on the scale of integrals 
over whole clouds? 

(iv) What is the intrinsic scatter in Ico as a given Ay? This cap¬ 
tures the degree to which this one dimensional approach represents 
a reasonable shorthand for the complex geometry of real clouds. 

A useful goal to enable the sort of calculation we describe above 
would be a library of relations that capture the realistic spread 
among this relation in the Milky Way and Magellanic clouds and 
allows for an understanding of how the key features such as the 
threshold, scatter at fixed Ay, saturation level, and dependence on 
environment. 

Our knowledge of the PDF requires similar refinement. The 
PDF of Milky Way clouds at low column remains substantially un¬ 
known (Lombardi et al. 2015), leading to uncertainties in the func¬ 
tional form of the gas distribution. Similarly, direct knowledge of 
the PDF of clouds in other galaxies is almost totally lacking. In 
the coming years both advances will help our understanding of the 
physics of CO emission substantially. 

Finally, a careful handling of the different phases of the ISM 
will improve our understanding of the situation. We have adopted 
a very direct observational approach in this paper, simply consider¬ 
ing all dust and CO in each « 10 pc Magellanic Cloud beam. Via 
comparison with Hi it should eventually be possible to model line 
of sight contamination by dust unassociated with the cloud, though 
this introduces subtleties regarding what dust is relevant for shield¬ 
ing. An improved analytic treatment of the H 1 -H 2 breakdown could 
also help refine the calculation. 


5 SUMMARY 

We show that at 10 pc resolution the relationship between dust col¬ 
umn expressed as visual extinction and CO intensity appears sim¬ 
ilar in the low metallicity Magellanic Clouds and the Milky Way. 
This agreement across a range of metallicity supports the theoreti¬ 
cally motivated view of shielding by dust as the dominant factor in 
determining the distribution of bright CO emission. To show this, 
we use surveys of CO emission from the Large (Wong et al. 2011) 
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and Small (Rubio et al., in preparation) Magellanic Cloud. We com¬ 
bine these with estimates of Av based on Herschel infrared maps 
from the HERITAGE survey (Meixner et al. 2010). We compare 
the Magellanic Cloud measurements to highly resolved Milky Way 
observations for two clouds, matched resolution measurements for 
local molecular clouds, and high latitude CO and dust emission as 
seen by Planck. 

Our measurements are consistent with an approximately uni¬ 
versal relationship between CO intensity and dust extinction within 
molecular clouds, though with only « 10 pc resolution we do not 
conclusively demonstrate such a relationship. Even for an approxi¬ 
mately universal relation, we still expect such a relationship to vary 
at second order due to changing geometry and environment. We 
show suggestive evidence for such a variation in the Large Mag¬ 
ellanic Cloud, where lines of sight with cooler dust temperatures 
show brighter CO emission at fixed Av- This could indicate that 
the weaker radiation field in these regions lowers the density of 
dissociating photons, allowing CO to emerge at fainter Ay. 

We discuss the implications of a nearly universal Ico-Ay re¬ 
lationship and suggest a simple, separable model for thinking about 
integrated CO emission from molecular clouds. In this picture, the 
PDF of column densities, the dust-to-gas ratio, the Ico-Av rela¬ 
tion, and the H 2 -H 1 boundary combine to determine the properties 
of a cloud but can be treated as separate problems. A number of 
studies have already considered parts of this problem as separable. 
Here we explore the implications for the CO-to-H 2 conversion fac¬ 
tor of a varying dust-to-gas ratio and fixed Ico-Av relation. Taking 
the PDF of local molecular clouds, we calculate the corresponding 
Av distribution for a range of dust-to-gas ratios and then predict 
the CO emission for each case. The result is a prediction for the 
variation of CO-to-H 2 conversion factor. 

Our empirically-motivated model predicts Xco oc Z~^ to 

in above about 1/3 Zq, in rough agreement with a variety of 
previous observational and theoretical studies. Our calculation also 
highlights the tenuous nature of CO as a tracer of molecular mass 
at metallicities even as high as 1I5Zq. At these metallicities, both 
the details of the H 2 -H 1 transition and the shape of the high end of 
the column density PDF will be extremely important to Xco ■ For 
a range of reasonable assumptions our calculations yield Xco that 
can scatter by as much as an order of magnitude at these metallici¬ 
ties. Future work will be useful to establish the functional form and 
variation of the Ico-Av relation and PDF within clouds, including 
those at low metallicity and in other Local Group galaxies. 


APPENDIX A: EFFECT OF BASELINE SUBTRACTION, A 
THICK GALACTIC HI DISK, AND 30-DORADUS ON THE 
7co-^y RELATIONS 

Al Baseline Subtraction 

In MAGMA DR2, a single linear baseline with magnitude of or¬ 
der a few mK has been subtracted from all spectra. In this paper, 
we carry out an additional additional zeroth order baseline subtrac¬ 
tion from the MAGMA data cube pixel-by-pixel. Our local base¬ 
line has mean magnitude « 1.2 mK with 39 mK rms scatter from 
position to position. Because we integrate the data cube over the 
whole velocity axis (266 channels with a channel width of ~ 0.5 
km s“^), this baseline will change the local integrated CO inten¬ 
sity by « 0.2 ± 5 K km s“^. In the left panel of Figure Al, we 
show the affect of our additional additional baseline subtraction on 
the Ico-Av relation in the LMC. The additional baseline correc¬ 
tion effectively removes the CO emission at low Av, which we 
thus interpret as likely artifacts. Overall, though, the relation does 
not change much and the baseline correction is almost irrelevant at 
high Av- 


A2 The Contribution of the Milky Way’s Thick Hi Disk to 
Av at High Galactic Latitudes 

We argued that the Ico/Av for the high latitude Milky Way is 
likely biased low because of contamination by dust associated with 
a long path length through the Milky Way Hi disk. To estimate the 
magnitude of this effect, we compute a simple correction factor by 
adopting a typical hydrogen nuclei number density n « 0.1 cm“® 
for the WNM component with a scale height of 400 pc at the lo¬ 
cation of the Sun (Table 1 in Kalberla 2003). Assuming a typical 
Galactic dust-to-gas ratio, the Av from this component would be 


. „ , _3 400 pc X 3.09 X 10^® cm pc ^ 

Al/= 0.1 cm ®- -A- - -1 - A- - 

sin(b) X 1.87 x 10^^ cm“^ mag~i 


(Al) 


where we convert column density to Av adopting a Galactic dust- 
to-gas ratio (X(H) = 5.8 x 10^^ E{B — V), Bohlin et al. 1978) 
and Rv = 3.1, as in Section 3.4. 

In the middle panel of Figure Al, we show Ico-Av relation 
in the high Galactic latitude lines of sight in the Milky Way after 
correcting for the 400 pc thick Milky Way Hi disk contribution to 
Av- Compared to Figure 6, the average trend of Milky Way high 
latitude lines of sight in Figure Al is much more closer to the aver¬ 
age LMC Ico-Av relation. 
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A3 30-Dor versus Non-30-Dor Ico-Av Relations 

In Section 3.3 we used the strength of interstellar radiation field 
traced by Tdust as a second parameter to test its effect on Ico-Av 
relation in the LMC. We also considered dividing the Large Magel¬ 
lanic Cloud into two regions, 30-Dor and Non-30-Dor (where 30- 
Dor region is defined as a rectangular box surrounding the Hll re¬ 
gion), to see if there is any systematic effect of star bursting envi¬ 
ronment on the Ico-Av relation. The Ico-Av relations in 30-Dor 
and Non-30-Dor are shown in the right panel of Figure AL Unlike 
the case for Tdust, we do not see any notable differences between 
the Ico-Av relations in 30-Dor and non-30-Dor, except at very 
low Av lines of sight. We suspect that the weird behaviour of the 
relation at low Ay for 30-Dor region is likely arising from unstable 
baseline corrections towards CO faint lines of sight. Considering 
the large dispersion of Ico at a given Av, we conclude that there 
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Figure Al. (lefty. Effect of baseline subtraction, before (red circles) and after (grey downward triangles; shifted slightly toward right in x-axis for easier 
comparison) our correction. Additional baseline subtraction effectively removes the fixed level of CO emission at low Ay while leaving the relation unaffected 
at high Ay. (middle)'. Ico-Ay relation in the high latitude Milky Way lines of sight after subtracting the contribution of 400 pc thick Hi disk to Ay (black 
triangles). The median Iqq at a given Ay in the high latitude Milky Way lines of sight is now closer to the average EMC value (red circles), compared to 
Figure 6. (right): Ico-Ay relation in the LMC near 30 Doradus (red upward triangles) and outside 30 Doradus (blue downward triangles). Although there 
are differences at low Ay, there is no obvious strong effect associated with being near 30 Doradus above our completeness limit at Ay Ri 0.8 mag. (A colour 
version of this figure is available in the online journal.) 


is no noticeable difference in the Ico-Ay relation between the two 
regions. 


APPENDIX B: STACKED SPECTRA FOR THE LMC 

Our working assumption in the main text of the paper is that simi¬ 
lar physics operate at small scales in molecular clouds in the LMC, 
SMC, and Milky Way. We interpret the coincidence of LMC, SMC, 
and Milky Way data at low resolution in Ico-Ay space to support 
this. A consequence of this conclusion is that the very low line- 
integrated intensities seen in the Magellanic Clouds are often a re¬ 
sult of dilution by our 10 pc beam. 

As a basic test of this, we stack the MAGMA spectra in bins 
of Ay and examine their peak temperature and line width in each 
bin. If beam dilution is the dominant physics in setting the observed 
stacked intensity, then we expect the peak temperature of the spec¬ 
trum to vary monotonically with Ay and integrated intensity while 
the line width shows no clear trend. In this case the diminishing in¬ 
tensity simply represents averaging similar spectra with more and 
more empty space within the beam. We show that this is the case 
in Figure Bl, which plots CO spectra from LMC in bins of Ay. 
Before averaging, the velocity axis of the CO line has been shuf¬ 
fled using Hi emission (Kim et al. 2003) as a template (following 
Schruba et al. 2011). 

In Figure Bl, we show the stacked spectrum in each Ay bin 
as grey filled histograms. We fit a single Gaussian to each stacked 
spectrum to derive the peak brightness temperature (Tpk) and full 
width at half maximum (FWHM) of the CO line. We overplot resid¬ 
uals from the Gaussian fit as red diagonally hatched histograms. 
Figure Bl clearly shows that the peak brightness temperature in¬ 
creases as a function of Ay except in the last Ay bin where a mis¬ 
match of CO and Hi velocity is apparent and leads to lower Tpk 
than the real value. On the other hand, the fit line width does not 
have a clear trend. It increases at low Ay and decreases at high Ay. 
Therefore, it appears that the peak main beam temperature drives 
the rise of integrated CO intensity as a function of Ay in the LMC. 

The peak values lie far below the expected kinetic tempera¬ 
tures of the molecular gas, so we expect that the dominant factor in 


Figure B1 is a changing level of beam dilution. Changing excita¬ 
tion may play some role as well; Pineda et al. (2008) found that the 
CO excitation temperature increases as a function of visual extinc¬ 
tion in the Perseus molecular cloud (See Figure 10 of their paper), 
ranging from 5 K at low Ay to 20 K at high Ay. Increased heating 
associated with molecular peaks might also play a role, either due 
to enhanced star formation activity (e.g., Heiderman et al. 2010) or 
efficient photoelectric heating (Hughes et al. 2010). 

Large beam dilution effects has been seen in the Magellanic 
Clouds as early as Rubio et al. (1993) and Lequeux et al. (1994), 
when it was invoked to explain the weak CO emission in the SMC. 
It has also been noted in CO observations of the LMC using SEST 
telescope (Garay et al. 2002; Kutner et al. 1997). We show that 
this effect remains strong even at 10 pc resolution. The direct re¬ 
sult of this calculation is that the Ico-Ay relation in the Magel¬ 
lanic Clouds must be carefully interpreted. It will not straightfor¬ 
wardly encode features like saturation and an Ay threshold (see 
also Wong et al. 2011). However, once data are matched in scale, 
we expect that its similarity to Milky Way results does indicate sim¬ 
ilar physics operating at sub resolution. 
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